The Photon Wavefunction: a covariant formulation and equivalence with QED 
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We discuss the limits of the photon wavefunction (PWF) formalism, which is experiencing a 
revival in these days from the new practical applications in photonics and quantum optics. We build 
a Dirac-like equation for the PWF written in a manifestly covariant form and show that, in presence 
of charged matter fields, it reproduces the standard formulation of (classical) Electrodinamics. This 
shows the inconsistency of the attempts to construct a quantum theory of interacting photons, based 
on the so called photon wavefunction approach, alternative to standard QED. PWF formalism can 
then be used to provide an easier description of the propagation of free photons, when the photon 
number remains fixed in time. 
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INTRODUCTION 

The problem of writing a wavefunction for the pho- 
ton takes its origins from the first attempts of quantiz- 
ing the electromagnetic field, since the birth of Quan- 
tum Mechanics. Because of the localization problem 
for the photon, the definition and even the existence 
of a wavefunction for the photon is still controversial, 
as it cannot always give a complete description of the 
system, like that provided by Schrodinger equation for 
non~relativistic massive particles. The first attempts 
can be found in the unpublished notes by Majorana []|, 
where the quantum states of the electromagnetic field 
were tentatively described by using the language of first 
quantization in the form of a Dirac-like equation, ob- 
tained from the Riemann-Silbcrstein (RS) formulation of 
Maxwell equations 2, 3, 4]. Dirac equation was formu- 
lated to describe the evolution of the relativistic electron, 
a particle with non-zero rest mass, ?i/2 spin, and elemen- 
tary charge e. Weyl equations instead describe massless 
neutral spinors. Finally Majorana extended the Dirac 
equation also to particles with arbitrary spin, in a more 
general infinite-spin component formalism [5|. This ap- 
proach for the quantization of the Electromagnetic field 
in the first quantization language is justified by the fact 
that Maxwell equations present an intrinsic mathemati- 
cal structure similar to that of a quantum wave function 
in relativistic theory and, conversely, the same procedure 
followed by Dirac to write the relativistic equation for the 
electron can be used to derive Maxwell equations. New 
recent experiments with single, double and many-photon 
sources and also with entangled states, where the photon 
number is small and remains fixed^enewed the interest 
in the wavefunction of the photon 0, S B S K El, [II • 
This revival of interest raised some, in our opinion too 
optimistic, hopes that this approach could open new per- 
spectives for alternative quantum descriptions of pho- 
tons, even not free [13[. The Majorana-like equation for 
the PWF can be considered as a consistent Quantum 
Mechanics of a free photon only, even if modulo the well 



known difficulty of its non localization \\4\ . For a system 
of relativistic particles in interaction (and in particular 
of interacting photons) Quantum Mechanics cannot be 
used since it implies action-at-distance forces, incompat- 
ible with Relativity. Quantum Field Theories are, in fact, 
mandatory. 

In this letter, after rewriting the Majorana-like equa- 
tion of the photon wavefunction in a manifest covari- 
ant form, we discuss the lagrangian that reproduces this 
equation, to be used as the starting point for getting a 
Quantum Field Theory for (interacting) photons. Then 
we show that this lagrangian is equivalent to that of clas- 
sical Electrodynamics so that after quantization it leads 
to standard QED. 

PWF AND COVARIANT FORMULATION 

Following Majorana formulation, without loosing in 
generality for a particular choice of helicity state, or 
of multiplicative constants, one defines the Riemann- 
Silbcrstein vector 



E 



±iB 



and Maxwell Equations in the vacuum become 



V-r = 0, iVxF = ±- 



- dt 



(1) 



(2) 



By using the correspondence principle, p <-> p = ^ihV 
{i — 1,2,3), F here represents the wavefunction of the 
photon, leading to the wave equation =F^^F-f ip x F = 
0, while V • F = is the transversal! ty of the fields with 
respect to the propagation direction, namely p • F = 0. 
By introducing the 3x3 complex matrices 



s^ = I -i 

i 

one obtains a Dirac-like equation 

d 




''dl^ 



HF, 



(3) 



with H = ±cs • p and s = {sx-,Sy,Sz) [li la III- This 
Hamihonian has eigenvalues ±cp, 0. The eigenvalue is 
forbidden by the transversality condition. 

There has been a debate in the literature on the in- 
terpretation of the negative energy state with eigenvalue 
—cp. The interpretation of the states with eigenvalues 
±cp as states of positive energy and helicity ±1 is dis- 
cussed in [6|. A simple way to understand this point is 
the following: since in our case the observables of Energy 
and Helicity commute, we can interpret the generator of 
translations in time, H in ^ (after diagonalization) as 
the product of the Hamiltonian H — cp with the helicity 
operator A. The wave equation now has energy eigenval- 
ues always positive and the negative eigenvalue of H is 
due to the negative helicity value. As for the trasversality 
condition V • F = 0, let us recall that it is at the origin 
of the non localization of the PWF. Indeed the Hilbert 
space where F{x) lives, is the space of modulo, square 
functions 4'{x) that satisfy the condition V • = and 
therefore, the position operator X(j){x) — x(f){x) is not an 
operator of this space, since x4){x) does not satisfy the 
transversality condition. 

The first step to a covariant formulation for F is the 
isophormism between the algebras of the group SL{2C), 
of unimodular 2x2 matrices in the complex field and the 
(proper orthocronus) Lorentz group, 5*0(1,3), of 4 x 4 
(pseudo-)orthogonal real matrices that leave invariant 
the Minkowsky metrics rj^'^ = diag{l, —1, —1, —1). Since 
the two algebras are isomorphic, the two groups satisfy 
a local isomorphism that, extended to a global one, be- 
comes a 2 ^ 1 homomorphism. Consider the matrix 



Xq + X3 Xl — iX2 
Xi + iX2 Xf) - X3 



(4) 



built with the space-time coordinates, X{),xi,X2tX^ and 
transform X by an 5L(2C) trasformation as x' = A~^xA, 
where A e 5L(2C), so that detx' = detx. Since detx — 
Xq — |xp = ry^'^x^Xj,, the transformation leaves invariant 
the four dimensional interval s^ = rj^'^Xi^Xi, and therefore 
induce a Lorentz transformation A on x^, but to both 
±A correspond the same Lorentz transformation (i.e. the 
homomorphism is 2 — » 1 ). The group SL(2C) has two 
inequivalent, fundamental representations called (jjO) 
and (O, i) that can be also considered as two spinorial 
representations of 50(1,3) (chiral and antichiral Wcyl 
spinors). 

The elements of the vector space on which these repre- 
sentations operate are (j)a{x), {a = 1,2) for the (^,0) 

and 6 = e""( 



{a, a — l,2)for (O,^), such that 



2 



(l)'^(x') = AJ<I)i3{x) and </> (x') = {A'' - ')"' ^f (x) where 
A,A'<-^ €SL{2C) 

Tensorial products of the fundamental representations 
of the Lorentz group give rise to higher dimensional rep- 
resentations, divided in two main classes: Tensorial rep- 
resentations that derive from the product of even times 



the fundamental representations (^, 0) and/or (0, ^), and 
Spinorial representations coming from the product of 
odd ones. The simplest tensorial representations are 

ijj^ = (j)a® X'^, 1paf3 = (j)a® X/3, V'd/j ='4'a® X.f3, with 

0a, X/3 G (57 O), 0a, X^ e (O, i). Written in its symmet- 
ric and antisymmetric parts, V'a/s = V'[q/3] + ^(ai3)j then, 
V''[q/3] = Aeq^, has the only degree of freedom A, where 
e"^ is the 2x2 completely antisymmetric tensor. The 
symmetric part has instead 3 independent components. 
A similar decomposition holds for ip^g- The constant ten- 
sors £q/3, e^g and their inverses e"^ , e°-^ can be used to 
rise and lower the spinorial indices. Let us define the two 
by two matrices {cr'^)^ and (ct'')q'' where a^ ^ cr'^ = 1, 
CT* = — CT* and CT* are the Pauli matrices. The tensorial 
representations ip^, ipap and ij:^^ can be expressed in 
terms of these matrices. 






(5) 
(6) 
(7) 



but the first of these equations, that defines a four- vector 
"0^. Similarly, 



1 



0'[mH =^[^.u\±i7^£^^l'P'y^• 



per 



(8) 



Both 0;+ and ^ have 3 independent components 
and ip+^ = 1!;+^, -0;7^ = -VV' where ipf,^° = 
j^e/iiypCT'0'"^ is the dual tensor of -0^,^. 

One can verify easly that ip~ and ip~^ do not contribute 
to 0'q,/3 and to ipaa respectively so that 






(9) 
(10) 



Therefore V'(q/3) describes the self-dual part, ipii,] and 
Tpi^m the antisclf-dual part,'0r , , of a six-components 

double-antisymmetric tensor -0(^1/] € (^ ^,0) © (O, ^ i). 
The Faraday electromagnetic tensor F^^ such that 
F°' = —-Ei and F^^ = —e^^'^Bk is a double antisym- 
metric tensor and its self-dual and anti self-dual parts 
I F'^^ , Fr, 1 can be written in covariant spinor notation 



- ('^['^TT'^l P.+ 



i^(o/3) = (^"^^'^'^[;.])a/3 



F 



(d/3) 



(ctI^ct-'If, 



[p,.]M/3 



(11) 
(12) 



Where F^^^^ e (i i, O) and F, .^) G (O, i i), but 



i^(.« « ((ct[V1)F+) ^^^(ct'^+)„,, (13) 



'a/3 ■ 



(./3) « ((-'V1)F„-t)^^^(ct*^-),^ (14) 



and 

Ft = ^± zB, (15) 

so that i^(Q/3) and F,^m are the positive and negative he- 
hcity wavefunctions of the photon in covariant notations. 
To write the Dirac-Hke equation for the photon wave 
function, consider 

{^''d^)jF^p^^ = {a^d^<j''a^)o.aF+^^0 (16) 

and its complex conjugate 



(^''9^)/^ 



(/3d) 



(a^a^aV^)„<iF-;, = 0. (17) 



When eqn. (fTH]) is saturated with (cr'^)"'^ one gets, for 
T = 0, the first equation in ([2|) and, for t — i, the second 
equation in ^ that is equation ^ . The same results are 
obtained saturating (fT7|) with (ct'^)"". 

Equation ([3]), toghether with the transversality condi- 
tion V • F = 0, is equivalent to free Maxwell equations. 
This leads to speculate that this approach could be taken 
as the starting point for a new quantum description of, 
even not free, photons. However when it is rewritten in 
covariant form, (|16p . it becomes completely clear that it 
describes just free Maxwell equations in a different nota- 
tion. That leaves little room to the speculations previ- 
ously mentioned. 

Let us add some further considerations to stress this 
point ever more. As noted in the introduction, a rela- 
tivistic quantum theory with interactions must be neces- 
sarly a local QFT. The recipe to write the (free) classical 
field lagrangian density, to be quantized, is to look at 
the classic action that yields the Schroedinger equation 
of the Quantum Mechanics of the single particle (even- 
tually supplemented with local interaction terms) and 
quantize this classical action according to the canonical 
rules. This procedure has been named, quite improperly, 
second quantization. 

For instance the lagrangian density of a free elec- 
tron, derived from the Dirac equation, is Ce = 
^ (17^9^ — m/4) ^, which is invariant under the global 



gauge transformation ^ 



^, where A is the global 



gauge parameter. To extend this trasformation to a lo- 
cal one with gauge parameter A (a;) one must introduce a 
gauge field A^ that trasforms a,s A^ ^ A 
lagrangian density becomes 



d^X and the 



^e = * (i-f^id^ - ieA^) - mh) *. 



(18) 



Now we have to search for a lagrangian density that 
gives rise to (fT6|) as its Eulero-Lagrange field equations. 
Since (|16|) transform covariantly as a four-vector one 
needs a four- vector, let say written in spinor notations, 
j^aa _ j^^f^^fi-^aa ^ j^ lagrangian density that reproduces 
(fTB)) by varying A is £ = a{Aa^df_i)°'^Faf3 where a is 



a normalization constant. Notice that this lagrangian is 
invariant under the gauge trasformation A^ — > A^ -|-9^A. 

If one adds to the lagrangian C the Dirac lagrangian Ce 
it is right to identify, as anticipated with the notations, 
the real four-vector in A with the A^ in (jlSp since with 
this identification the field equations for >1^ yield cor- 
rectly the current term, j'^ = e'^j^'^ in the right hand 
side of the first group of Maxwell equations (or equiva- 
lently a non linear term j* in the r.h.s. of ^ and a non 
linear term j° in the r.h.s. of V • F = 0). However the 
field equations obtained varying F(q,^) are (^j^A^])"*" = 
that imply A^^ = modulo a gauge transformation. But 
C has a serious drawback: it is not real. 

By adding to C the complex conjugate counterpart, 
the action becomes 



/ = / -[{Aa'^d^rf'Fe.p + {Aa^^d^r^F^^] + / C, 



= I aiA^d^Ff^n + J ^e 



(19) 



where A = A\a . Unfortunatly (1191) reproduces only 
the first group of Maxwell equations d^^F^^, = ej„ (with 
a = 3/4). 

A possible cure of this desease could be to allows A 
and A to become complex that is ^4 = {A\ + iB\)a^, 
A = {Ax - iBxp^ {Afj, and B^ real). Now dUD (with 
current) are reproduced correctly but there is a unac- 
ceptable doubling of degree of freedom. Indeed now the 
action ITOl) becomes 



/ = JiA^d^F'^'' + B^e^'^P'^d^Fp, + Ce) (20) 

The second group of Maxwell equations e^^'''"^duFp„ = 
implies F^^ = ^{d^A^ — 9^^^). When this algebraic 
equation is used to remove F^^, in (PO]) the B^ field drops 
out and the action becomes / = J{d[pA^]d^f^A'^^ + Ce) 
Even worse, if one defines A^^'^ = ^{A^ ± A^), then 

1= /(9[^4|^9[^A''i'+^-9[X7^a[''A''i(-)+/:^) (21) 



so that, after quantization, one of the two "photons" de- 
scribed by the gauge fields A^^ has negative metric and 
therefore the action ((2T|) is inconsistent. 

The only consistent way to cure these problems is to 
add to the lagrangian density in (fT9l) the term 



3 
32' 



^0 = -i^^iF'^'^Fp^ + F"^V^„) = -F^^F^., (22) 



to obtain 



1 = 1 Ctot = - J aM A^lF^, - j ^F'^'^Fp, (23) 
+ / ^(^7^(9^ - ieAp) ~ mh) ^ 



The field equation for i^^^ identifies F^j^ with 9[^Ay]: 
-F)ii/ = -l/2{d^A^ ~ d^A/^,) = d[f_iA^]. Since this is an 
algebraic equation it can be used to replace F^^ in the 
lagrangian with the following result, 



/ = 



^[^A.jai'^A'^i 



+ / *(i7^(9^-ieA^)-m/4)* 



(24) 



eqn. ([24|l is the standard action of classical electrody- 
namics and after quantization it gives rise standard QED. 
In the presence of charged matter fields, this formulation 
reproduces only the standard formulation of (classical) 
Electrodinamics. This clearly shows the inconsistency 
of the attempts to construct a quantum theory of in- 
teracting photons, based on the so called photon wave- 
function approach, alternative to standard QED. There 
is a perfect correspondence between PWF and QED only 
when photons are free, non interacting and when the pho- 
ton number remains constant during the evolution of the 
field, with the problems of the photon localization. 



DISCUSSION AND CONCLUSIONS 

In this note we have presented and discussed the equiv- 
alence of the PWF formalism with that of standard 
Quantum Electrodynamics. PWF can only describe sce- 
narios where the photons are free, non interacting and 
maintain a constant number during their evolution. No 
absorption and/or emission of photons can be directly de- 
scribed by the formalism of PWF. The equivalence is set 
by a manifestly covariant version of the so called PWF 
equation. Moreover, on the basis of this covariant formu- 
lation we have motivated the statement that the photon 
wave function approach, at the second quantization level, 
cannot give anything else than the standard Quantum 
ElectroDynamics . 

The correspondence set between PWF and QED for 
free photons is useful to shed some light in the Orbital 
Angular Momentum (QAM) of the photon and the PWF 
formalism [H, [H, [13, [III . 0AM of light is deeply con- 
nected with the vorticity of the E-M field and with the 
creation of optical vortices. Recently the RS vector, 
which is the basis with which the PWF is built, was used 
to describe the E-M field vorticity [H, S, HJ [111 . Rie- 
mann Silberstein vortices are defined by 



F(r,i) ■F{r,t) = 0. 



(25) 



The loci of points satisfying this condition are lines in 
space, the phase of the field is singular surrounded by 
zones where the phase gradient vector is circulating. 
Laguerre-Gaussian beams are particular cases in which 
the field has spatial symmetry and the RS vortex lines are 
stationary. Anyway exact solutions of electromagnetic 



waves carrying angular momentum have been recently 
described in by using the momentum representation and 
then were cast in terms of PWF in the RS formalism 
|22| that, by using the correspondence here discussed, is 
simply equivalent to the description obtained with QED 
by quantizing the field in paraxial approximation that, 
at the single photon level, represents the probability am- 
plitude of finding a photon in a certain eigenstate of mo- 
mentum, helicity and QAM |23l|. t hat is clearly not an 
intrinsic property of the photon[24J. The equivalence of 
the two formulations can be easily set by expressing the 
PWF in terms of the photon annihilation and creation 
operators, using the vector potential A in the RS vector, 
but this goes beyond the purpose of this work. 
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By definition, the intrinsic properties of a particle are 
those that do not depend on the choice of a reference 
frame, i.e the rest mass, the electric charge and spin. 



